Saclay T02/038 
LBNL-50102 
heFph /0204142| 



Fine Structure Constant Variation 
from a Late Phase Transition 



Z. Chacko"'", C. Grojean"^ and M. Perelstein" 

" Theory Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
^ Department of Physics, University of California, Berkeley, CA 94720, USA 
^ Service de Physique Theorique, CEA-Saclay, F-91191 Gif-sur-Yvette, France 
zchackoSthsrv . Ibl . gov , gro j eanOspht . saclay . cea . f r , meperelsteinSlbl . gov 

Abstract 

Recent experimental data indicates that the fine structure constant a may be 
varying on cosmological time scales. We consider the possibility that such a vari- 
ation could be induced by a second order phase transition which occurs at late 
times (z ~ 1 - 3) and involves a change in the vacuum expectation value (vev) 
of a scalar with milli-eV mass. Such light scalars are natural in super symmetric 
theories with low SUSY breaking scale. If the vev of this scalar contributes to 
masses of electrically charged fields, the low-energy value of a changes during the 
phase transition. The observational predictions of this scenario include isotope- 
dependent deviations from Newtonian gravity at sub-millimeter distances, and (if 
the phase transition is a sharp event on cosmological time scales) the presence of a 
well-defined step-like feature in the a{z) plot. The relation between the fractional 
changes in a and the QCD confinement scale is highly model dependent, and even 
in grand unified theories the change in a does not need to be accompanied by a 
large shift in nucleon masses. 



A recent analysis of quasar absorption spectra at redshift z ~ 0.5-3.5 indicates 
that the fine structure constant a is changing with time: 

anew - apa.t ^ 5a ^ ^ ^^g^ (J) 

a now a 

While theorists have considered the possibility that the fundamental constants are time- 
dependent for a long time (starting with Dirac [2]), it is not clear how the result fits 
into the current field-theoretic picture of elementary particle physics. In the Standard 
Model (SM), all the coupling constants run with energy. Since the temperature of the 
Universe changes with time, this implies that the effective values of the couplings are 
also changing. This effect, however, cannot be used to explain the data, since a does 
not run at energies below the electron mass me ~ 500 keV, corresponding to redshifts 
of order 10^. Thus, it appears that the time variation of a reported in pP cannot be 
explained without appealing to physics beyond the Standard Model. 

A priori, string theory should be well suited to accommodate this data. Indeed, 
it is well known that in string theory any coupling constant is promoted to a vacuum 
expectation value (vev) of a scalar field such as the dilaton or some other modulus. If 
this scalar field is extremely light, m ~ 10~^^ eV, its expectation value could be still 
evolving in the recent past (or even today.) However, it is not clear how such a low 
scalar mass can be stabilized against radiative corrections. Moreover, the superlight 
scalar will mediate a long-range, isotope-dependent force. Such forces are subject to 
severe constraints from accurate tests of the equivalence principle (for a recent review, 
see 0) and other fifth force experiments Finally, it has been pointed out |7j that 

the energy density associated with the rolling scalar field is large enough to overdose the 
Universe if the parameters of the model are chosen to accommodate the result ([Q). (This 
problem seems to be tightly connected with the usual cosmological constant problem, 
and may become a non-issue once that problem is resolved [8..) 

In this paper we propose an alternative explanation of the time variation of the fine 
structure constant that is potentially experimentally distinguishable from the slowly 
rolling dilaton scenario. Our framework is based on the observation that the experi- 
ments measure the value of the fine structure constant at the very low energy scale 
(about 10 eV) set by atomic physics. This infrared value of the constant, which we will 
denote by a = a(0), is related to its ultraviolet value a(A) by the renormalization group 
equation: 

a ayh) 2tt ~^ rrii 

where mo < mi < m2 < . . . < mjv are the masses of all electrically charged particles, 
commonly referred to as "thresholds" (mg = rrie, m-i = m^, etc.) and bi is the one- 
loop beta function between the scale mj_i and rrii. For notational convenience, we have 
identified mAr+i with the cutoff of the theory A. With our conventions, 6 < for an 
asymptotically free theory. It is clear from that the low-energy value of a can change 
even if its fundamental, short- distance value a(A) remains fixed, provided that some of 
the thresholds move. It is this possibility that we would like to investigate in this paper. 



1 



A change in a induced by a small variation of thresholds is easily obtained from 0: 



-=-f (3) 

In QED, bi — = —qfrii < 0, where is the electric charge of the particle that 
decouples at rrii, and is its degeneracy. To induce the change of a at 10~^ level, we 
need 

(4) 

i * 

Notice that the masses have to increase (mnow — '^past > 0) in order to reproduce the 
correct sign in the evolution observed by Ref . . The requirement (0} is not in obvious 
contradiction with other observations. While the bounds on the variation of electron 
and proton masses at 2; ~ 2 — 3 are rather tight j^, there are no constraints on the 
masses of other charged particles {e.g. muon or r) at these redshifts. 

What physical mechanism could lead to a change in the mass of a charged particle? 
In many theories masses of charged fermions are determined by vacuum expectation 
values of electrically neutral scalar fields. (For example, lepton and quark masses in the 
Standard Model are proportional to the Higgs vev.) In the expanding Universe, scalar 
vevs can be time-dependent. In particular, we will consider models in which in the early 
Universe the vev of a certain scalar S vanishes due to thermal effects. As the Universe 
cools down, a phase transition occurs, during which the scalar acquires a vev. As a 
result, the masses of electrically charged states coupled to S will change, leading to a 
change in the low-energy value of a according to Q • 

Let us determine the features of the zero-temperature scalar field potential V{S) 
which are necessary to explain the data 1 . We assume that in the early Universe 
5 = 0. It is easy to see that V'{0) has to vanish; however, 5 = could be either a 
minimum or a (local) maximum of V{S). In the latter case, the vacuum can still be 
stable, provided that a thermal mass term is generated by interactions of the S field 
with the surrounding plasma [TUj : 

Kherm ^ T"^ S'^ . (5) 

For the moment, we will simply assume that such a mass term is generated; we will later 
comment on the conditions under which this is the case. As the temperature drops to 
its critical value, characterized by 

Tl - -V"{Q), (6) 

the vacuum at S* = becomes unstable, and the system undergoes a phase transition. 
In this transition, the field changes its vev until it reaches the nearest minimum of 
V{S), which we will denote by ^i. We require that the phase transition occurs at 
2; ~ 1 — 3, when the temperature of the Universe T is of order 10~^ eV. Eq. ^ then 
implies that the mass parameter of the S field in the zero-temperature lagrangian is of 
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the same order, 10~^ eV. In a generic non-supersymmetric theory, such a low mass scale 
is unstable against radiative corrections, and extreme fine tuning would be required to 
explain it. In supersymmetric theories, however, this scale can be radiatively stable, 
provided that the supersymmetry breaking scale is at or below about 10 TeV and 
the breaking is only communicated to S via Planck suppressed contact operators^. Such 
a low value of \/F is phenomenologically viable if the breaking is communicated to 
the visible sector fields (that is, the Standard Model fields and their superpartners) by 
gauge mediation Thus, we have identified a large class of models in which a phase 
transition could occur at low redshifts, without any fine-tuning. 

The phase transition will not result in a change of the fine structure constant unless S 
is coupled to at least one electrically charged field. On the other hand, radiative stability 
of the 5" mass parameter would be destroyed by any direct renormalizable coupling of S 
to the Standard Model fields or their superpartners. To avoid this problem, we introduce 
an additional vectorlike chiral superfield Q, which is charged under U{l)em (and possibly 
other SM gauge groups.) This field is coupled to S via 

C = j d?e{M + yS)QQ, (7) 

where M is a mass for the Q and Q fields, and ?/ ~ 1 is a Yukawa coupling. It is easy to 
see that three-loop diagrams involving Q and Q renormalize the S mass parameter by 
an amount 

These corrections do not destroy radiative stability provided that M ^ 10^^ GeV. Note 
that Eqs. (jU and ((Tj) imply that Si/M ~ 10~^ is necessary to accommodate the data, 
so the vev of the 5* field after the transition is required to be very large. Si ^ 10^^ 
GeV. The huge hierarchy between the vev of the field and its mass is radiatively stable 
because of supersymmetry and implies that S" is a modulus. 

Let us return to the question of whether a thermal mass term of the form (jSj) is 
actually generated at high temperatures. Since the couplings of 5* to the visible sector 
fields are non-renormalizable, they cannot generate such a term. Moreover, since is a 
modulus, its self-interactions are extremely weak: the quartic coupling can be estimated 
as A ~ rn^/Sf ~ 10~^^. The thermal mass term generated by these interactions is 
proportional to A and is therefore negligible. To generate a thermal mass of the right 
order of magnitude, S needs to have substantial couplings to some other light fields 
Y. These couplings will not destabilize the 10"'^ eV mass scale, provided that the 
fields Y, like S itself, are only sensitive to supersymmetry violating effects through 
Planck suppressed operators. (Together with S, these fields form a "hidden sector" of 
the theory.) Although the fields of the hidden sector decouple from the visible sector 
fields very early on, they can maintain thermal equilibrium among themselves, at a 

^The fact that the energy scale corresponding to the current critical density of the Universe could 
arise naturally in supersymmetric theories was emphasized in Ref. jllj . 
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temperature which is identical (up to order one factors related to the multiplicity of 
states at decoupling) to the visible sector temperature. The interactions of S with this 
"hidden" thermal plasma are responsible for generating the mass term ©. 

The natural value of the energy density difference before and after the phase transi- 
tion is given by 



and given the above constraint we obtain AV > eV^. Naively, such a huge energy 
density seems to lead to two severe problems, which would make our scenario incom- 
patible with standard cosmology. First, one could argue that the cosmological constant 
cannot be tuned away both before and after the phase transition, and rapid inflationary 
expansion should occur in at least one of these two periods. Second, during the phase 
transition the field S is expected to undergo coherent oscillations, slowly decaying due to 
Hubble expansion. With our parameters, the energy density in these oscillations is large 
enough to overdose the Universe. Requiring that our scenario without any additional 
ingredients be consistent with standard cosmology leads to a bound on the variation of 
a similar to the one found in Ref. 0: 6a/a ^ 10~^^, well below the values reported 



A simple way to restore the consistency of our scenario with standard cosmology is 
to simply fine-tune the shape of the potential. Such fine tuning can be used to get rid of 
the large contribution to the vacuum energy before the phase transition, Eq. 0. Once 
this is done, the energy of coherent oscillations is automatically low enough to avoid 
overclosure. The amount of fine tuning involoved is similar to what is required in the 
conventional models with a rolling dilaton field [?]. 

It is tempting, however, to speculate that the difficulties of our scenario could be 
avoided without any fine tuning once the cosmological constant problem is resolved. 
Indeed, from the effective field theory point of view, it is reasonable to expect that the 
mechanism that solves the cosmological constant problem will guarantee the absence of 
inflation regardless of the (constant) vacuum energy density^. During a phase transition, 
such a mechanism would adjust itself to cancel the new value of the energy density. Any 
mechanism which has this feature would resolve the flrst of the two problems confronting 
our scenario. Furthermore, the energy in the coherent oscillations of the modulus fleld 
is related to the difference in the vacuum energies before and after the transition. It 
is the large value of this difference that leads to the second problem of our scenario. 
If the large difference in vacuum energies is cancelled by the adjustment mechanism, 
conservation of energy implies that the amount of energy in the coherent oscillations 
will be sufficiently small to avoid overclosure. Of course, we emphasize that no explicit 
example of an adjustment mechanism for the cosmological constant problem is currently 
known, and therefore the above discussion is necessarily rather speculative. If such a 
mechanism is proposed in the future, the question of whether our scenario is viabile 
without fine tuning should be reexamined within a more concrete framework. 

^This point of view was advocated, for example, in Ref. 




(9) 



in p. 
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In the context of a grand unified theory a change in a would seem to require a 
variation of the QCD confinement scale. This, in turn, would result in a shift of the 
hadron masses. In fact, if the observed variation of a is attributed to the change of 
the short- distance unified couphng acuT, the fractional change in the proton mass can 
be predicted [HIISI (see also (TH] for an estimate in the change of the deuteron binding 
energy) and turns out to be about 40 times larger than the fractional change in a. This 
result is in mild contradiction with the bound obtained from measuring the value of 
= mf./mp [9 in the same range of redshifts, 2; ~ 2 — 3. In our framework, however, it 
is straightforward to evade this prediction. This possibility is perfectly compatible with 
grand unification. For example, consider a supersymmetric SO(IO) grand unified theory 
in which doublet-triplet splitting is achieved using an adjoint field S with block diagonal 
vacuum expectation value (0, 0, 0, r2, where rs are the Pauli matrices [inillZl- This 
pattern of vevs breaks SO(IO) down to SU(4) x SU(2)iX U(1)/3r, where the U(l) is 
the diagonal generator of SU(2)jj. The superpotential of the theory necessarily contains 
a Planck suppressed operator of the form ST,QQ/Mpianck, where Q and Q are in the 
fundamental of SO(IO), and S' is a singlet. Because of the pattern of vevs of the field 
S, the singlet 5* acquires Yukawa couplings of the form (|7j) to the doublets in Q and Q 
but not the triplets. When S goes through a phase transition, the low energy values of 
the 311(2)2, and U(l)y gauge couplings are affected, but the SU(3) gauge coupling (and 
therefore the QCD confinement scale) remains unaffected at one loop order, and the 
resulting change in the proton mass will therefore be suppressed. However, there will be 
changes in both the electron mass and the proton mass arising from radiative corrections 
to these quantities involving electromagnetic and weak interactions. Moreover, since su- 
persymmetry breaking is mediated to the visible sector fields through gauge interactions, 
the Higgs mass and therefore its vacuum expectation value will be altered by the change 
in the values of the SU(2) and U(l) gauge couplings. Both these effects will in turn 
give rise to a change in /i of order Sfi/fi ~ 6a/a. Although this is below the current 
experimental bound, it may be possible to detect this effect in future experiments. 

The observational predictions of our model depend on the time T it takes to com- 
plete the phase transition. Since the dynamics of the transition is necessarily strongly 
influenced by the unknown adjustment mechanism for cosmological constant, it is not 
possible to estimate T reliably. Let us consider two cases. In the first case, T is smaller 
than the age of the Universe at the time when the transition began, so that the transi- 
tion is a sharp event on cosmological time scales. Then, the low-energy value of the fine 
structure constant should be time-independent for z ^ 1. This prediction is consistent 
with the null results from the laboratory search ^H] for the time variation of a, as well 
as the geological bound ^H] from the Oklo natural nuclear reactor. Moreover, in this 
case we would expect improved measurements of a{z) to show a sharp feature (a "step") 
in the z ~ 1 — 3 region. The other possibility is that the time T is so large that the 
transition is still not completed today. In this case, it is much harder to distinguish our 
scenario from the more conventional picture of a slowly rolling dilaton field. 

Apart from the quasar absorption spectra measurements, the only currently available 
bound on the variation of a on cosmological time scales comes from nucleosynthesis |2Uj : 



5 



our model easily satisfies this constraint (|Aa|/a < 10^^ — 10^"^). The expected preci- 
sion [21] in the cosmic microwave background measurements (|A«|/« < 10^^ — 10^^) 
will not be enough to rule out our scenario. 

An observational consequence of the rolling dilaton mechanism of changing a is the 
existence of a new, isotope-dependent long-range force jHEllSl- In our scenario, the mod- 
ulus 5* (in the true vacuum) has a mass of about 10~^ eV and mediates a Yukawa force 
with a range of order 0.1 mm. Since there are no renormalizable couplings between 5* and 
any of the Standard Model fields, this force is of (approximately) gravitational strength. 
Currently, the strongest bounds on such forces come from tabletop precision tests of the 
Newton's law at submillimeter distances These experiments require that the range 
of the extra force be less than about 0.2 - 0.3 mm. Clearly, this bound can be satisfied 
in our model without significant fine tuning. If our scenario is indeed realized in nature, 
improved precision gravitational tests should observe isotope-dependent deviations from 
the Newton's law in the near future. 

Above, we have assumed that S* = is a local maximum of the zero-temperature 
potential V[S)^ and the phase transition is second order. It is also possible that S* = 
is a local minimum of the potential. If there is another minimum with lower energy, 
at 5 = 5*1, a first order phase transition could occur. During this phase transition, 
bubbles of the true vacuum nucleate and start growing rapidly. Eventually, the bubbles 
coalesce and the transition is complete. A classic calculation (neglecting the effects of 
the expansion of the Universe) of the lifetime of a false vacuum [SB] yields 



where A is the natural scale of the potential V{S\ ISV = V{0) — V{Si) is the energy 
splitting between the true and false vacua, and 



The phenomenological constraints on this scenario are somewhat different from the case 
of a second order transition discussed above. According to (fTUj). the lifetime of the false 
vacuum can be much larger than A without major fine tuning. Thus, it is possible 
for a late phase transition to occur even if all the scales in V^(5') are of order TeV or 
higher. This in turn allows direct couplings of S to the (electrically charged) visible 
sector fields. Note, however, that the absence of rapid infiationary expansion before 
(or after) the phase transition would naively require AV^ < (10^^ eV)^. This bound 
would not only necessitate fine tuning of the energy splitting itself, but also lead to 
unacceptable large values of r (much larger than the age of the Universe today) unless 
the entire shape of V{S) is finely tuned. On the other hand, just like in the case of 
second order phase transition discussed earlier, this bound could be nullified once the 
physics responsible for solving the cosmological constant is taken into account. 

To summarize, we have suggested that the apparent time variation of the low energy 
fine structure constant reported in is a result of a late phase transition involving 




(10) 




(11) 



6 



changing the vacuum expectation value of a hght modulus. We have shown that such 
a late phase transition can occur in a large class of supersymmetric models without 
any fine tuning. This is in marked contrast with a more conventional rolling dilaton 
mechanism of changing a, which requires superlight (Hubble- scale) moduli whose mass 
is generally not radiatively stable even in the absence of gravity. Just like in the rolling 
dilaton picture, producing a change in a large enough to explain the data PP in our 
scenario is challenging, and naively seems incompatible with cosmological observations. 
In both pictures, this problem could be avoided by fine tuning the potential of the scalar 
field. We have speculated that our scenario could be compatible with standard cosmol- 
ogy without fine tuning, provided that the cosmological constant problem is resolved by 
an adjustment mechanism. Clearly, this question would have to be reexamined if an ex- 
plicit adjustment mechanism is proposed in the future. While we are not able to reliably 
describe the dynamics of the phase transition, we emphasize that the transition could be 
a sharp event on cosmological time scales. In this case, improved measurements of the 
redshift dependence of a will discover a well-defined step-like feature, distinguishing this 
scenario from the rolling dilaton picture. Another prediction is the isotope-dependent 
deviation from Newtonian gravity at sub-millimeter scales. Finally, we have found that 
in our framework the relation between the change in a and the corresponding change 
in the QCD confinement scale is highly model dependent. In particular, even in su- 
persymmetric grand unified theories it is straightforward to construct models which are 
consistent with current bounds on the variation of the ratio of the electron mass to the 
proton mass at 2; ~ 2 — 3. 

Acknowledgments 

We thank C. Csaki, E. Dudas, N. Kaloper, M. A. Luty, H. Murayama, A. E. Nelson 
and J. P. Uzan for fruitful discussions. M.P. wishes to acknowledge the hospitality of 
Service de Physique Theorique, CEA-Saclay, where this work was initiated. Z.C and 
M.P. are supported by the Director, Office of Science, Office of High Energy and Nuclear 
Physics, of the U. S. Department of Energy under Contract DE-AC03-76SF00098. 

References 

[1] J. K. Webb, M. T. Murphy, V. V. Flambaum, V. A. Dzuba, J. D. Bar- 
row, C. W. Churchill, J. X. Prochaska and A. M. Wolfe, Phys. Rev. Lett. 
87, 091301 (2001) iastro-ph700T 2539|; see also J. K. Webb, V. V. Flambaum, 
C. W. Churchill, M. J. Drinkwater and J. D. Barrow, Phys. Rev. Lett. 82, 884 
(1999) astro-ph/9803165|; V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Phys. 
Rev. Lett. 82, 888 (1999). 

[2] P. A. Dirac, Nature 139, 323 (1937); for an historial perspective, see F. Dyson, "The 
fundamental constants and their time variation", in Aspects of Quantum Theory, 



7 



eds A. Salam and E. Wigner. 



[3] T. Damour, |gr-qc/0109063| 

[4] J. D. Bekenstein, Phys. Rev. D 25, 1527 (1982). 

[5] G. R. D vali a nd M. Zaldarriaga, Phys. Rev. Lett. 88, 091303 (2002) 
|hep-ph/0108217| . 

[6] K. A. Olive and M. Pospelov, Phys. Rev. D 65, 085044 (2002) | hep-ph/0110377| . 

[7] T. Banks, M. Dine and M. R. Douglas, Phys. Rev. Lett. 88, 131301 (2002) 
ihep-ph/0112059j . 

[8] P. Langacker, G. Segre and M. J. Strassler, Phys. Lett. B 528, 121 (2002) 
lhep-ph /0112233( . 

[9] A. Y. Potekhin, A. V. Ivanchik, D. A. Varshalovich, K. M. Lanzetta, J. A. Bald- 
win, G. M. Williger and R. F. Carswell, Astrophys. J. 505, 523 (1998) 
|astro-ph/9804116| . 

[10] S. Weinberg, Phys. Rev. D 9, 3357 (1974); L. Dolan and R. Jackiw, Phys. Rev. D 
9, 3320 (1974). 

[11] N. Arkani-Hamed, L. J. Hall, C. Kolda and H. Murayama, Phys. Rev. Lett. 85, 
4434 (2000) |astro- ph/000"5TTl|. 

[12] M. Dine, W. Fischler and M. Srednicki, Nucl. Phys. B 189, 575 (1981); S. Dimopou- 
los and S. Raby, Nucl. Phys. B 192, 353 (1981); L. Alvarez-Gaume, M. Claudson 
and M. B. Wise, Nucl. Phys. B 207, 96 (1982); M. Dine and A. E. Nelson, Phys. 
Rev. D 48, 1277 (1993) fhi^ph/9303230 ; M. Dine, A. E. Nelson and Y. Shirman, 
Phys. Rev. D 51, 1362 (1995) |hep-ph/9408384 ; M. Dine, A. E. Nelson, Y. Nir and 
Y. Shirman, Phys. Rev. D 53, 2658 (1996) hep-ph/9507378 |. 

[13] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper and R. Sundrum, Phys. Lett. B 480, 
193 (2000) hep-th/0001197|; S. Kachru, M. Schulz and E. Silverstein, Phys. Rev. 
D 62, 045021 (2000) |hep-th/0 001206 ; C. Csaki, J. Erlich and C. Grojean, Nucl. 
Phys. B 604, 312 (200i nhep-th/00l2T43| . 

[14] X. Calmet and H. Fritzsch, |hep-ph/0112Tl0l H. Fritzsch, |hep-ph/020lT98 

[15] T. Dent and M. Fairbairn, Nucl. Phys. B 653, 256 (2003) |hep-ph/0112"279| . 

[16] G. R. Dvah and S. Pokorski, Phys. Lett. B 379, 126 (1996) |hep-ph/9601358|. 

[17] Z. Chacko and R. N. Mohapatra, Phys. Rev. Lett. 82, 2836 (1999) 
[hep-ph/9810315j . 



8 



[18] J. P. Turneaure, C. M. Will, B. F. Farrell, E. M. Mattison and R. F. C. Vessot, 
Phys. Rev. D 27, 1705 (1983); J. D. Prestage, R. L. Tjoelker and L. Maleki, Pliys. 
Rev. Lett. 74, 3511 (1995) 

[19] A. I. Shlyakhter, Nature 264, 340 (1976); T. Damour and F. Dyson, Nucl. Phys. 
B 480, 37 (1996) liep-ph/9606486|; Y. Fujii, A. Iwamoto, T. Fukahori, T. Ohnuki, 
M. Nakagawa, H. Hidaka, Y. Oura, P. M511er, Nucl. Phys. B 573, 377 (2000) 
fhep -ph/9809549j . 

[20] E. W. Kolb, M. J. Perry and T. P. Walker, Phys. Rev. D 33, 869 (1986); B. A. Camp- 
bell and K. A. Olive, Phys. Lett. B 345, 429 (1995) hep-ph/9411272|; L. Bergstrom, 
S. Iguri and H. Rubinstein, Phys. Rev. D 60, 045005 (1999) astro-p h/9902157j ; 
P. P. Avelino et al. Phys. Rev. D 64, 103505 (2001) |astro-ph/0102144 |. 

[21] S. Hannestad, Phys. Rev. D 60, 023515 (1999) |astro-ph/9810102l; M. Kaplinghat, 
R. J. Scherrer and M. S. Turner, Phys. Rev. D 60, 023516 (1999) astro-ph/9810133|. 
P. P. Avelino, C. J. Martins, G. Rocha and P. Viana, Phys. Rev. D 62, 123508 (2000) 
|aitro-p h/0008446 . 

[22] C. D. Hoyle, U. Schmidt, B. R. Heckel, E. G. Adelberger, J. H. Gu ndlach, D. J. Kap - 
ner and H. E. Swanson, Phys. Rev. Lett. 86, 14 18 (200 1) |h ep-ph/ 00 11014 . 
E. G. Adelberger [EOT- WASH Group Collaboration] , |hi^ x/0202008| 

[23] S. Coleman, Aspects of symmetry. (Cambridge Univ. Press, 1985.) 



9 



